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Cn ' Abstract. This paper provides a brief review of the history of our understanding and knowledge of black holes. Starting 

;_( ' with early speculations on "dark stars" I discuss the Schwarzschild "black hole" solution to Einstein's field equations and 

^1^, the development of its interpretation from "physically meaningless" to describing the perhaps most exotic and yet "most 

■^^ . perfect" macroscopic object in the universe. I describe different astrophysical black hole populations and discuss some of 
their observational evidence. Finally I close by speculating about future observations of black holes with the new generation 

CO ' of gravitational wave detectors. 



> THEORETICAL CONSIDERATIONS 

^. 

^D ■ This paper summarizes a talk presented at the Albert Einstein Century International Conference, held in Paris, France, 

in July 2005 to mark the centennial of Einstein's "miracle" year 1905. Strictly speaking, black holes are a consequence 
of another Einstein miracle that did not occur in 1905 and still had to wait ten years, namely general relativity. With 

^^ ' his theory of special relativity, however, Einstein had laid the ground work already in 1905 for his theory of general 

C^ . relativity. And even more strictly speaking, speculations on black holes predate even the special theory of relativity by 

"^ ' over a century. 



P ! Early Speculations 



j^ , In the late 1700's John Mitchell [1] in England and Jean Simon Laplace [2] in France independently realized that 

celestial bodies that are both small and massive may become invisible' . The basis for this speculation is the observation 
that the escape speed 



^, (1) 

where M and R are the stellar mass and radius, is independent of the mass of the test particle. Within Newton's particle 
theory of light it seems quite reasonable that this should also apply to light, in which case light can no longer escape 
the star if the escape speed exceeds the speed of light, 

Vesc > C. (2) 

Evidently this happens when 

c'^R 
GM>—, (3) 

meaning that stars with a large enough mass and a small enough radius become "dark". Laplace went on to speculate 
that such objects may not only exist, but even in as great a number as the visible stars. With the demise of the particle 
theory of light, however, these speculations also lost popularity, and dark stars remained obscure until well after the 
development of general relativity. 



An excellent review of the histoiy of our understanding of compact objects can be found in [3]. 



General Relativity and the Schwarzschild solution 

In 1915 Albert Einstein published his field equations of general relativity [4], 

Gat = ^TlGTah- (4) 

One could argue that superficially this equation may not be all that different from the Newtonian field equation 

v2^=47rGp. (5) 

The left-hand-side of the Newtonian field equation (5) features a second derivative of the Newtonian potential ^, 
and the right-hand-side contains matter densities p. Quite similarly the Einstein tensor Gab on the left-hand-side of 
Einstein's field equations (4) contains second derivatives of the fundamental object of general relativity, the spacetime 
metric gab- To complete the analogy, the stress-energy tensor Tab on the right-hand-side contains matter sources. For 
the vacuum solutions in which we are interested in this article, the stress-energy tensor vanishes, T^t = 0. 

Unfortunately the Einstein tensor Gab contains many lower-order, non-linear terms, making Einstein's equations a 
complicated set of ten coupled, quasi-linear equations for the ten independent components of the spacetime metric gab- 
Clearly, it is very difficult to find meaningful exact solutions. Karl Schwarzschild, returning fatally wounded from the 
battle fields of World War 1, was nevertheless able to derive a fully non-linear solution in spherical symmetry within 
a year of Einstein's original publication [5]. Written as a line element ds^, the spacetime metric gab describing this 
solution is 

ds^ ^gabdx''dx'' = -i\- — — \dt^+i\- — — \ dr^ + r^de^ + r^sm^ed(p^. (6) 

This solution is the direct analog of the Newtonian point-mass solution 

*^2^ (7) 

r 

and describes the strength of the gravitational fields, expressed by the spacetime metric gab, created by a point mass M 
at a distance r. The relativistic Schwarzschild solution (6) is significantly more mysterious than its Newtonian analog, 
though. Particularly puzzling is the "Schwarzschild radius" 

2GM 

rs = —r- (8) 

at which the metric (6) becomes singular The existence of this singularity obscured the physical interpretation of this 
solution, and in fact Schwarzschild himself died believing that his metric was physically irrelevant. 

Other aspects also contributed to the fact that the astrophysical significance of the Schwarzschild solution (6) 
remained unappreciated for decades. For one thing it was not clear how such an object could possibly form, even 
though in 1939 Oppenheimer and Snyder [6] published a remarkable analytic calculation describing the "continued 
gravitational contraction" of a dust ball that leaves behind the Schwarzschild metric (6). This calculation, serving as 
a crude model of stellar collapse, made several simplifying assumptions: that the matter has zero pressure (so that it 
can be described as dust), that the angular momentum is zero, and that the spacetime is spherically symmetric. Critics 
maintained that in any realistic astrophysical situation none of these assumptions would hold, and that any deviation 
from these idealizations could easily halt the collapse, preventing continued gravitational contraction. The absence 
of angular momentum seems particularly troubling. Since the Schwarzschild solution does not carry any angular 
momentum it was completely unclear how any astrophysical object, which would necessarily carry some angular 
momentum, could collapse and form such a solution. 

Finally, no astronomical observations had revealed phenomena requiring a gravitationally collapsed object as an 
explanation. Given that at least some of the astronomical community had been quite reluctant to accept the much less 
exotic White Dwarfs as an astrophysical reality, it is not surprising that a solution to Einstein's field equations that 
became singular at finite radius was not immediately embraced as a celestial object. 

All three of these factors - a lack of understanding of the singularity at the Schwarzschild radius, general skepticism 
concerning gravitational collapse, and the absence of astronomical observations of gravitationally collapsed objects - 
resulted in the fact that the great astrophysical significance of Schwarzschild's solution (6) remained unappreciated for 
almost 50 years. All of this changed in the 1960's, the "golden age of black hole physics". 
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FIGURE 1. Cartoon of a spacetime singularity surrounded by an event liorizon at rj = 2M. Tiie liglit cones of an observer falling 
into the black hole become increasingly tilted, until they completely tip over at the event horizon. 

THE GOLDEN AGE OF BLACK HOLE PHYSICS 

The golden age of black hole physics was ushered in when advances in our theoretical understanding of the 
Schwarzschild geometry and gravitational collapse in general coincided with new astronomical observations of highly 
energetic objects that clearly pointed to a gravitationally collapsed object as their central engine. 

The Schwarzschild geometry 

To begin with, it became clear that the apparent singularity at the Schwarzschild radius rs (8) is simply a harmless 
coordinate singularity, not completely unlike the poles of a sphere when described in terms of longitude and latitude. 
This was demonstrated by Kruskal [7], who introduced a new coordinate system^ that remains perfectly regular at rs. 
Instead of being a mysterious singularity the Schwarzschild radius now emerged as a point of no return: a one-way 
membrane called "event horizon" through which nothing, not even light can leave the collapsed region inside. Inside 
this event horizon lurks a true singularity at which the curvature of spacetime becomes infinite^. 

Consider the unfortunate observer in the spacetime cartoon of Figure 1. As long as he is far away from the event 
horizon an "outgoing" light ray emitted away from the horizon can propagate toward larger distance almost as easily 
as an "ingoing" ray emitted toward the horizon can propagate toward smaller distance. At this point the observer's 
local light cone, which is the section of spacetime on which light propagates, is almost upright, and he can easily send 
signals to a buddy even further away. As our first observer approaches the event horizon, however, his light cones 
become increasingly tilted, and it becomes increasingly difficult for outgoing light rays to actually move outward. At 



^ Interesting, at least one other regular coordinate system had already been found by Painleve [8] and Gullstrand [9]. Lemaitre [10] similarly 

concluded that the geometry at the Schwai'zschild radius must be regular, but apparently these conclusions remained unnoticed. 

^ In a classical description of general relativity, that is. A self-consistent quantum theory of gravity will presumably resolve this singularity. 



the event horizon the Ught cone tips over. An "outgoing" Hght ray emitted exactly on the event horizon will hover at 
the Schwarzschild radius forever, and inside this point all light rays, even "outgoing" ones, immediately move inward. 
Nothing can emerge from the event horizon and everything inside will soon reach the singularity that lurks at the 
black hole's center. Our observer will cross the event horizon and reach the singularity in finite proper time. As seen 
by his buddy far away, his signals slowly fade away as he approaches the event horizon. Note also that the in-falling 
observer would not necessarily observe anything special occurring at the event horizon. Ultimately the observer will 
be torn apart by the increasingly strong tidal forces, but that may occur inside or outside the event horizon, depending 
on the black hole's mass. It is only through careful experiments with flashlights, studying the properties of outgoing 
light-rays, that the observer would be able to tell that he has crossed an event horizon. 

It is interesting to note that, by sheer coincidence, the location of the event horizon r^ (8) as expressed in 
Schwarzschild coordinates coincides exactly with the radius of a "dark star" (3) as determined by Mitchell and Laplace. 



Gravitational Collapse 

Almost simultaneously with the improved understanding of the Schwarzschild geometry, it became clear that grav- 
itational collapse is much more generic than was believed earlier. In 1963 Roy Kerr [11] discovered a generalization 
of the Schwarzschild solution (6) that carries angular momentum. With this discovery the argument that objects car- 
rying angular momentum cannot collapse gravitationally immediately collapsed itself. The first numerical relativity 
simulations also demonstrated that pressure may not be able to prevent gravitational collapse [12]. 

Further progress in our understanding of gravitational collapse arrived with a number of theorems. Perhaps most 
importantly Roger Penrose [13] showed that the formation of a spacetime singularity is generic after a so-called trapped 
surface has formed. Another set of theorems are collectively called the "no-hair" theorems^; in essence these theorems 
state that black holes have no distinguishing features. That is to say that Kerr's solution describing a rotating black 
hole is the only solution describing a rotating black hole - all stationary black holes are Kerr black holes, parametrized 
only by their mass M and the angular momentum J^. Black holes can be perturbed, but any perturbation is quickly 
radiated away, leaving behind a Kerr black hole. 

This is a truly remarkable statement: it means that the structure of black holes is uniquely determined by their mass 
and angular momentum alone, completely independently of what formed the black hole at first place. This observation 
led Chandrasekhar to eloquently conclude in his Nobel Lecture 

This is the only instance we have of an exact description of a macroscopic object. . . They are, thus, almost 
by definition, the most perfect macroscopic objects there are in the universe. 

Observational Evidence 

The technical advances during World War II lead to the development of radio astronomy in the post-war years. 
Within a few years several discrete sources had been detected, but with few exceptions the origin of these sources 
remained a mystery. It was generally believed that the sources were otherwise dark "radio stars" in our galaxy, since 
at extra-galactic distances they would have to be enormously energetic. This opinion started to change when the 
positioning of these "quasars" improved, and some were identified with galaxies. 

The true breakthrough arrived when lunar occultations of the radio source 3C273 lead to its identification with 
another object whose optical spectrum showed a redshift of z = 0.158 - clearly establishing it as an extra-galactic 
object. This realization is documented in the remarkable March 16 issue of Nature [17, 18, 19, 20]. The same volume 
contains a theoretical paper by Hoyle and Fowler [21], in which the authors conclude^ 

Our present opinion is that only through the contraction of a mass of 10^ — lO^Mg to the relativity limit 
can the energies of the strongest sources be obtained. 



** See, e.g., Israel [14], Carter [15], Robinson [16]. 

' Strictly speaking black holes can also carry a charge Q, but from an astrophysical perspective this is irrelevant since any charge would very be 
neutralized very quickly. 
See also Ginzburg [22]. 



These events ushered in "Relativistic Astrophysics" as a new field. As a sign of the time the first Texas Symposium 
on Relativistic Astrophysics was convened in December of 1963. The new field also attracted many new people into 
the field, including John Wheeler. In fact, it was John Wheeler who in 1967 coined the term black hole, marking the 
transition from speculative ideas on dark stars to the astrophysical reality of black holes^. 



ASTROPHYSICAL BLACK HOLES 

Clearly, we do not have any absolutely water-tight proof that black holes exist in our universe. However, we do have 
some extremely convincing evidence that makes black holes by far the most conservative explanation of the observed 
phenomena. 

Observations clearly point to two different populations of black holes. One of these populations are "stellar-mass 
black holes", which have masses in the order of 10 Mg; another group are "supermassive black holes", which have 
masses in the order of 10^ — IQ^Mq. I will discuss these two groups in more detail below. There is also some 
observational evidence for "intermediate mass black holes" of about 10^ — lO^M©, and finally there are speculations 
on "primordial black holes" that would be left over from the big bang. Even though these different kinds of black 
holes have vastly different masses, they are, from a mathematical perspective, the exact same kind of animal: a Kerr 
solution to Einstein's field equations with certain values for their mass and angular momentum. From an astrophysical 
perspective they differ not only in their mass and angular momentum, but also in how we can observe them and in how 
they form, i.e. in their evolutionary history. 



Stellar-mass black holes: Cygnus X-1 

The prime example of a stellar-mass black hole is Cygnus X-1**, which was first discovered in the data of the X-ray 
satellite Uhuru, launched on Dec 12 1970. The case for Cygnus X-1 as a black hole can be summarized as follows^: To 
begin with, Cygnus X-1 shows very short time variations in the X-ray signal (in the order of Af ^ 10 ms and less). This 
implies that Cyg X-1 is a very small object, in the order of /? < cA/ < 10^ m. It is also the unseen binary companion to 
a 9th magnitude supergiant star called HDE226868. The Doppler curve of this star shows that the binary has an orbital 
period of about 5.6 days; from the amplitude of the Dopplershift we can determine the mass function / to be about 
O.25M0 [25]. Combining this with the mass of supergiant stars we can derive a lower limit on the mass of Cyg X-1, 

Mcygx-iS3.3M0. (9) 

An independent argument involving some constraints on the binary's orbit arrives at a similar limit. 

Given its size, Cyg X-1 has to be a "compact object": a white dwarf, a neutron star, or a black hole. Even under very 
conservative assumptions both white dwarfs and neutron stars have maximum masses safely below the lower limit of 
Cyg X-l's mass. This leaves a black hole as the most likely explanation. 

Cyg X-1 is an example of a stellar-mass black hole, which form the end-point of the evolutionary cycle for 
massive stars. When massive stars run out of nuclear fuel they can no longer support themselves against gravitational 
contraction. For stars with masses larger than about 20 M© nothing can halt the subsequent gravitational collapse, 
which therefore leads to prompt black hole formation. For stars with masses between about 8 and 20 M© the collapse 
can be halted when the density of the compressed stellar material reaches nuclear densities. The resulting shock-wave 
launches a "core-collapse" supernova and leaves behind a newly formed neutron star. This neutron star may either 
remain stable, or it may form a black hole at a later time. A "delayed collapse" to a black hole can be triggered by 
a variety of mechanisms, including fall-back of matter and phase transitions in the neutron star interior. Finally, two 
neutron stars may collide. If the remnant exceeds the maximum allowed mass for rotating neutron stars this coalescence 
will also lead to the formation of a stellar-mass black hole^°. 



' I again refer to Israel [3] for a much more complete account of the developments described in this Section. 

See Oda [23] for a review and references. 

See Section 13.5 and Box 13.1 in Shapiro and Teukolsky [24] for a detailed discussion. 

Even if the remnant exceeds the maximum allowed mass for uniformly rotating neutron stars it may be stabilized temporarily by virtue of 
differential rotation [26]. Dissipation of the differential rotation, for example by magnetic coupling, then triggers a delayed collapse 
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FIGURE 2. Orbit of the star S2 around SgrA* (Figure from [28].). 



We currently know of about 20 confirmed black hole binaries'', but presumably that is only a tiny fraction of the 
total number of stellar-mass black holes in our own galaxy. 

Supermassive black holes: Sagittarius A* 

Perhaps the most convincing evidence for a black hole comes from the center of our own galaxy, Sagittarius A* 
[28, 29] . Observations of our galactic center in the near infrared reveal several stars that orbit a central object in bound 
orbits'^. Particularly compelling is the orbit of the star S2, about two-thirds of which has now been mapped with 
increasingly accurate positioning (see Figure 2). The emerging orbit is a Kepler orbit with a period P of about 15 years 
and a semi-major axis a of about 4.62 mpc. From Kepler's third law we can conclude that the enclosed mass is 



2^3 



M = 



4n a 



'.Ax KfMr. 



(10) 



S2's orbit also has a significant eccentricity of e = 0.87. At pericenter its distance to the central object is only 124 
AU. This implies that the central object harbors an enormous mass in a very small volume. The most conservative 
explanation for such an object - which also must have remained stable over the lifetime of the galaxy - is again a 
black hole. 

Sagittarius A* is an example of a supermassive black hole. There is some very convincing evidence for supermassive 
black holes in other galaxies as well - for example the maser observations from NGC 4258 (also known as M106) 
[30]. In fact, there is evidence that supermassive black holes lurk at the cores of most galaxies [31]. 

It is less clear exactly how supermassive black holes form. Many different routes may lead to the formation of 
massive black holes in active galactic nuclei'^, but which of these routes nature tends to take is still under debate. A 
constraint comes from the recent observation of quasars at redshift z w 6 in the Sloan Digital Sky Survey [33]. If these 
quasars are indeed powered by supermassive black holes, this implies that the latter must have formed very quickly in 



" See Table 4.1 in McClintock and Remillard [27]. 

'^ A beautiful movie animation of these orbits can be found at www.mpe .mpg. de/ir/GC. 

'^ See, e.g., the flow chart in Fig. 1 of [32]. 




FIGURE 3. The binding energy E/mc^ as a function of separation rc^ jGM of a test mass m in orbit about a Newtonian point 
mass M (left panel) and a Schwarzschild black hole of mass M (right panel). The solid lines denote contours of constant orbital 
angular momentum L. Extrema of these contours identify circular orbits, marked by the dashed lines. The circular orbits are stable 
if the extremum is a minimum, otherwise they are unstable. 



the early universe. One model that may account for that is accretion onto seed black holes that form in the collapse of 
first-generation (Pop. Ill) stars [34, 35]. 



Probing black hole properties 



Most observational evidence for black holes to date is based on the argument "a lot of mass in a tiny volume", which 
leaves black holes as the most conservative explanation. Also, most observations to date only provide information 
about the black hole mass, and not about the angular momentum. Clearly it would be desirable to go beyond that. 

A number of efforts are under way to find observational evidence for certain black hole characteristics. One such 
characteristic is the absence of a stellar surface. X-rays emitted by black hole candidates indeed shows some differences 
from that emitted by neutron stars. These differences can be explained in terms of an accretion flow hitting the stellar 
surface in the case of neutron stars, but freely falling through the horizon in case of the black holes''*. Another effort 
aims at resolving the event horizon of Sgr A*, which would show up as a "black hole shadow" since it absorbs all 
radiation emitted behind it. 

An interesting idea for the measurement of the black hole angular momentum J is based on the concept of the 
Innermost Stable Circular Orbit, or ISCO for short. Consider a test mass m in orbit about a point mass M. A 
circular orbit can be found by identifying an extremum of the test mass's binding energy E at constant orbital angular 
momentum L; a minimum corresponds to a stable circular orbit, while a maximum corresponds to an unstable orbit. 

In Newtonian physics the binding energy as a function of separation r is given by 

L2 



G M 



2r2' 



(11) 



where L ~ L/mc^. To find an extremum of the binding energy we differentiate with respect to r at constant L and set 
the result to zero. 



L2 



GMr 



(12) 



'"* See Najayan [36] for a review and references. 



which we recognize as Kepler's third law (since L = Q.r^ /c^). Evidently we can find a circular orbit for arbitrary r, 
and the equilibrium energy of these circular orbits results from inserting (12) back into (11), 

1 GMm 

EhmA = 7: • (13) 

2 r 

We can also verify that all these extrema are minima, meaning that the orbits are stable. Alternatively, we can identify 
the orbits graphically as in the left panel in Fig. 3. 

For a test mass in orbit about a Schwarzschild black hole the binding energy is given by 

mc^ W c^ r J \ f^ J J 

where r is the Schwarzschild radius. Taking derivatives we would find that circular orbits exist only for radii r > 
3GM/c^. Moreover, these orbits are stable only outside 

6GM 

nsco = — n-, (15) 

which therefore marks the ISCO for a test particle orbiting a Schwarzschild black hole (see also the right panel in 
Fig. 3). 

The presence of an ISCO is of great relevance for black hole observations because most of the emitted radiation is 
believed to originate from an accretion disk. In this accretion disk particles follow almost circular orbits as they spiral 
toward the black hole. Since circular orbits are unstable inside the ISCO, the accretion disk can only exist outside 
nsco- Accordingly, the Doppler shift of spectral lines is limited by the speed of matter at the ISCO. The key point is 
that the location of the ISCO depends on the black hole's angular momentum J; it is at risco ~ 6GM/c^ only for a 
Schwarzschild black hole with J — 0, and can get much closer to the event horizon for a spinning Kerr black hole. For 
a spinning black hole the accretion disk may therefore extend closer to event horizon, and the correspondingly higher 
speeds result in a greater broadening of the emitted spectral lines. Some results based on this idea have been reported 
in [37]. 

Even if successful, however, these techniques can only determine the global parameters M and J. Clearly it would 
be desirable to map out the local properties of the spacetime geometry around a black hole. Our best chance of doing 
that comes with gravitational wave observations. 

GRAVITATIONAL RADIATION 

Sources of gravitational radiation 

Maxwell's equations predict that accelerating charges emit electromagnetic radiation. For an electric dipole, the 
emitted power is given by the Lamor formula 

2 .. .. 

^em = —^didi, (16) 

where di is the electric dipole moment, di its second time derivative, and where we sum over repeated indices. 

Einstein's equations similarly predict that accelerating masses emit gravitational radiation. One might expect that 
the power emitted from a gravitational wave source is given by an expression similar to Lamor's formula. However, 
the analog of the electric dipole moment is the mass dipole moment, and its first time derivative the total momentum. 
For an isolated system the total momentum is constant, so that the second derivative of the dipole moment vanishes - 
there is no gravitational dipole radiation. The first - and often dominant - contribution to gravitational radiation comes 
from the quadrupole term. The equivalent of the Lamor formula in general relativity is therefore 

iow = ^ ^ i^jk^jk) (17) 

Here J^ij is the reduced quadrupole moment 
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FIGURE 4. Cartoon describing tiie inspiral of a black hole binary. 

the bracket () denotes averaging over several characteristic periods of the system, and the triple dot denotes the third 
time derivative. For a strong signal of gravitational radiation we evidently need large and rapidly changing quadrupole 
moments, which brings to mind binary systems. 

To estimate Lqw for ^ binary system we evaluate J^jf; for a Newtonian binary and find 

32 C* mV^ 
Lgw = ^^—t-. (19) 

Here M = mi +m2 is the binary's total mass, jj. = m\m2/M the reduced mass and a the semi-major axis, and we have 
also used Kepler's law to eliminate the orbital frequency H. Inserting numbers we find that the only binary systems 
that emit appreciable amounts of gravitational radiation have huge masses and very small binary separations a. The 
most promising candidates are therefore compact binaries consisting of black holes or neutron stars. 

The loss of energy due to the emission of gravitational radiation leads a shrinking of the binary's orbit. To see this 
we compute 

da dE/dt _ Lqw _ 6A ( G M^l^^^l^\ 



dt dE/da dE\)ind/da 

where 



r.2 



(20) 



£.,.^i^ (21) 

2 a 

is the generalization of (13). Computing the loss of angular momentum we would also find that that the emission of 
gravitational radiation leads to a reduction of the binary's eccentricity, i.e. to a circularization of its orbit. 

A binary inspiral then proceeds as illustrated in Fig. 4. Presumably the binary starts out at a large binary separation. 
The emission of gravitational radiation leads to a continuous decrease in the binary separation, and also to a decrease in 
the binary's eccentricity. At sufficiently late times we may therefore approximate the binary orbit as circular, except for 
the slow inspiral. As the binary separation shrinks, both the amplitude and the frequency of the emitted gravitational 
wave signal increases. This leads to the typical "chirp" signal sketched in Fig. 4. In analogy to our discussion of 
point masses orbiting a single black hole, the binary will at some point reach an ISCO. Inside this separation it is 
energetically favorable for the binary companions to abandon circular orbits, and instead plunge toward each other 



and merge. In the final "ring-down" phase the remnant will settle down quickly into an axisymmetric equilibrium 
object. 

Even for the most promising sources of gravitational radiation any astrophysical signal that we might hope for is 
going to be extremely weak. To see this we estimate the effect of a gravitational wave on the spacetime metric g^h- 
If the spacetime is almost flat, we may write gab = "Hab + hab, where r\ab is the flat Minkowski metric and where hah 
is the small perturbation that we observe as a gravitational wave. For a perturbation caused by a distant gravitational 
wave source we have 

Here the symbol TT denotes the "transverse traceless" part of the corresponding tensors, s is the distance from the 
observer to the gravitational wave source, and the quadrupole moment J^,- is evaluated at the retarded time t — s/c. We 
again evaluate this term for a Newtonian binary, and estimate the gravitational wave amplitude to be in the order of 

\ G , , G^ M M rs rs 
s c^ c^ s a s a 

where we have used Kepler's third law to eliminate ii and where rs is the Schwarzschild radius (8) of the mass M. 
Evidently we expect the strongest gravitational wave signal htj when the binary's semi-major axis a is small, i.e. close 
to the ISCO where a « rs (if we neglect factors of a few). 

Consider, for example, a stellar-mass compact binary that is coalescing somewhere in the Virgo cluster. For such a 
binary, rs is in the order of a few km, and the distance s is about 15 Mpc or a few 10^" km. We then have h « 10^^**. 
Since the corrections htj are added to the background metric, which is of order unity, the relative change in the 
spacetime metric is exceedingly small, and the proposal to measure these perturbations remarkably ambitious. 

The weakness of any gravitational wave signal that we may realistically expect has two important consequences 
for the prospect of observing them: we need extremely sensitive gravitational wave detectors, and we need accurate 
models of any potential sources to aid in the identification of any signal in the noisy output of the detector. We will 
discuss both of these aspects in the following two Sections. 

Gravitational wave detectors 

A new generation of gravitational wave observatories is now operational. Basically, these instruments are gigantic 
Michelson-Morley interferometers. The two LIGO observatories in the US, for example, have an arm-length of 4 km. 
Somewhat smaller instruments exist in Italy (VIRGO), Germany (GEO), and Japan (TAMA), and another observatory 
is currently under construction in Australia (ACIGA). The idea is that a gravitational wave htj passing through these 
interferometers will slightly distort the relative length of the two perpendicular arms. Tracking these distortions with 
the help of laser interferences as a function of time should reveal the passing gravitational wave signal. 

The ground-based detectors mentioned above all have arm-lengths in the order of a kilometer, which makes them 
sensitive to the gravitational radiation emitted from stellar-mass black holes'^. A space-based gravitational wave 
antenna LISA with an arm-length of several million kilometer is being planned; this instrument would be sensitive 
to gravitational radiation emitted from supermassive black holes. 

As we have seen in the previous section an astrophysical gravitational wave signal that we might realistically expect 
will lead to a tiny perturbation of the spacetime metric and hence to only a minuscule distortion in the arm-lengths 
of the interferometers - in fact, only a tiny fraction of the size of the nucleus of a hydrogen atom. It is therefore 
extremely difficult to reduce the dominant sources of noise - seismic, thermal and photon shot-noise - and increase 
the signal-to-noise ratio so that an astrophysical source can be identified unambiguously. 

It is therefore truly remarkable that the LIGO collaboration has recently achieved its design goal of a strain 
sensitivity exceeding 10^^^ for a large part of its frequency range (see Fig. 5). Given this sensitivity, and assuming a 
signal-to-noise ratio of eight, LIGO is now able to detect binary neutron stars to a distance somewhat greater than 10 
Mpc. For binary black holes with a slightly larger mass this range is also slightly larger, and already includes the Virgo 
cluster. 

Whether or not the current LIGO observatory, LIGO I, will be able to detect a binary black hole system - or any other 
source of gravitational radiation - depends on how often such binaries coalesce in our or our neighboring galaxies. 



In geometric units a solar mass equals 1 .48 km. 
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Strain Sensitivities for the LIGO Interferometers 
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FIGURE 5. Strain sensitivities for the LIGO Hanford (Washington) interferometer (Figure courtesy the LIGO Scientific Collab- 
oration). 

Estimates for binary neutron star systems can be based on the statistics of known binaries'^, while estimates for 
binary black hole systems, which have never been observed, are typically based on population synthesis calculations '^. 
According to these estimates it is not completely impossible that LIGO I will observe a compact binary, and it seems 
almost certain that an advanced LIGO II observatory, as it is currently being planned, will see many such sources. 

We finally point out that a gravitational wave detector tracks the amplitude of a gravitational wave at one particular 
point in space, but, unlike a telescope, does not produce an image. A single detector therefore cannot position any 
source. Two detectors can position a source to a ring in the sky, and with the world-wide network of detectors there is 
hope that any potential source can be positioned to within a reasonable accuracy. 



Source modeling 

Different approximations can be used to model the inspiral of a compact binary in its different phases. For the 
initial inspiral, while the binary separation is sufficiently large and the effects of relativistic and tidal interactions 
sufficiently small, post-Newtonian point-mass calculations provide excellent approximations'**. The very late stage 
can be modeled very accurately with the help of perturbation theory'"^. Neither one of these approximations can be 
used in the intermediate regime around the ISCO during which the binary emits the strongest gravitational wave signal. 
The most promising tool for the modeling of this dynamical phase of the binary coalescence and merger is numerical 
relativity^". 



'« See, e.g. [38, 39]. 

'^ For example [40]. 

'* See, e.g., [41], as well as L. Blanchet's article in this volume. 

" E.g. [42]. 

^^ See [43] for a recent review. 



Numerical relativity calculations typically adopt a 3+1 decomposition. In such a decomposition the four- 
dimensional space M is carved up into a foliation of three-dimensional spatial slices E, each one of which corresponds 
to an instant of constant coordinate time t. Einstein equations for the four-dimensional spacetime metric gah can then 
be rewritten as a set of three-dimensional equations for the three-dimensional, spatial metric Yij on the spatial slices 
E, as well as its time derivative. 

The general coordinate freedom of general relativity has an important consequence for the structure of these three- 
dimensional equations. Since there are three space coordinates and one time coordinate, we can choose four of the 
ten independent components of the symmetric spacetime metric gab freely. But that means that the ten equations 
in Einstein's equations (4) cannot be independent - otherwise the metric would be over-determined. Four of the ten 
equations in Einstein's equations must be redundant. In the framework of a 3+ 1 decomposition these four equations are 
constraint equations that constrain the gravitational fields within each spatial slice E, while the remaining six evolution 
equations govern the time-evolution of Yij from one slice E to the next. The equations are compatible, meaning that 
fields that satisfy the constraint equations at one instant of time continue to satisfy the constraints at all later times if 
the fields are evolved with the evolution equations. This structure of the equations is very similar to that of Maxwell's 
equations, where the "div" equations constrain the electric and magnetic field at any instant of time, while the "curl" 
equations govern the dynamical evolution of the fields. 

Finding a numerical solution to Einstein's equations usually proceeds in two steps. In the first step the constraint 
equations are solved to construct initial data, describing the gravitational fields together with any matter or other 
sources at some initial time, and in the second step these data are evolved forward in time by solving the evolution 
equations. 

One of the challenges in constructing initial data results from the fact that the constraint equations determine only 
some of the gravitational fields. The remaining fields are related to the degrees of freedom associated with gravitational 
waves, which depend on the past history and cannot be determined from Einstein's equations at only one instant of 
time. Instead, these "background fields" are freely specifiable and have to be chosen before the constraint equations 
can be solved. The challenge, then, lies in making appropriate choices that reflect the astrophysical scenario one wishes 
to model. 

For the modeling of binary black holes we would like to construct initial data that model a binary at a reasonably 
close binary separation outside the ISCO as it emerges from the inspiral from a much larger separation. We expect 
such a binary to be in "quasi-equilibrium", meaning that the orbit is circular (except for the slow inspiral) and that the 
individual black holes are in equilibrium. A number of different approaches have been pursued, but currently the best 
approximations to quasi-equilibrium black holes are the models of Cook and Pfeiffer [44]. Probably there are ways 
to improve these data, in particular as far as the background fields and the resulting "gravitational wave content" are 
concerned^', but even without these improvements the current models are probably excellent approximations. 

While initial data describing binary black holes have been reasonably well understood for a number of years, 
progress on dynamical simulations of black holes has been much slower. Until recently, these simulations had been 
plagued by numerical instabilities that made the codes crash long before the black holes had done anything interesting. 
The past year, however, has seen a dramatic break-through in this field, and by now several groups can perform reliable 
simulations of the binary black hole coalescence, merger and ring-down. 

The first announcement of such a calculation came from Pretorius [48]. His approach differs from the more 
traditional numerical relativity calculations in that he integrates the four-dimensional Einstein equations (4) directly, 
instead of casting them into a 3+1 form^^. The left panel of Fig. 6 shows a gravitational wave form from one of his 
recent simulations, starting with the initial data of [44]. At early times some noise is seen, which is either related to 
the imperfect initial data themselves or the imperfect matching of the numerical methods used in the initial data and 
the evolution. The noise disappears quickly, and leaves behind a very clean waveform, tracking the inspiral through 
several orbits, through the plunge and merger, until the newly formed remnant settles down into a single Kerr black 
hole. 

Shortly after Pretorius' announcement a number of other groups announced similarly successful calculations 
[51, 52, 53]. All of these later simulations do adopt a 3+1 decomposition of Einstein's equations, and cast these 
equations in the BSSN form [54, 55]. The right panel of Fig. 6 shows the trajectory of the black holes in the calculation 
of Campanelli et al. [47]. Campanelli et al. [51, 47] and Baker et al. [52] adopt a particularly simple method that treats 



^' One approach is to match the background fields to the post-Newtonian approximation of the prior inspiral [45]; another promising suggestion is 

a "wave-less" approximation proposed by Shibata et al. [46]. 

^^ The equations are expressed in a way, however, that leads to particularly desirable mathematical properties; see [49, 50]. 
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FIGURE 6. The left panel shows the gravitational wave form Re(v/4) from a binary black hole coalescence as recently computed 
by Pretorius (private communication). The right panel shows the trajectory of the black hole horizons in the calculation of 
Campanelli et al. [47]. 

the black hole singularities as "punctures" and avoids having to excise the black hole interior from the numerical grid 
[56, 57]. In the meantime several follow-ups on these calculations have already been published, including simulations 
of binaries with mass-ratios different from unity [58, 59] and spinning black holes [60]. 

Especially comparing with the situation just a year ago, it is truly remarkable and reassuring that different groups 
using independent techniques and implementations can now carry out reliable simulations of binary black hole 
coalescence and merger. It may still be a while until results from these simulations can be used to assemble a catalog 
of realistic wave templates to be used in the analysis of data from gravitational wave detectors, but the past year has 
certainly seen a huge step forward in that direction. 



SUMMARY 



Black holes may well be the most fascinating consequence of Einstein's theory of general relativity. Similarly 
fascinating is the development of our understanding of black holes. Speculations on so-called "dark stars" actually 
predate both special and general relativity by over a century, but it nevertheless took almost half a century after 
the pubUcation of general relativity and Schwarzschild's derivation of his famous solution until black holes became 
generally accepted as an astrophysical reality. 

By now we have very convincing observational evidence for both stellar-mass and supermassive black holes. So 
far these observations only constrain the black hole mass, and clearly it would be desirable to probe the local black 
hole geometry in addition to the global parameters. There is hope that we will be able to do that in the near future 
with the new generation of gravitational wave detectors. The LIGO observatories have recently achieved their design 
sensitivities, enabling us to detect inspiraling binary black holes to a distance of approximately the Virgo Cluster. The 
next generation Advanced LIGO will improve the sensitivity by over a factor of ten, which increases the event rate 
by over a factor of thousand. With the recent advances in numerical relativity we are also much closer to producing 
theoretical gravitational wave templates, which will aid both in the identification of gravitational wave signals and in 
their interpretation. There is hope, then, that we can discuss detailed gravitational wave observations of binary black 
holes by the time we celebrate the centennial of Einstein's general relativity. 
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